function f = ind_shares(expdelta,expmu1,QI,ind_market)

ns=size(expmu1,2);
expon=repmat(expdelta,1,ns).*expmu1; 
denomi=(expon'*QI)'+1;      
denom=denomi(ind_market,:);    

f=expon./denom;
